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SUMMARY 

As  part  of  a  NAVAIRDEVCEN  Independent  Research  (IR)  program  on  the  aerodynamic 
characteristics  of  lifting  membranes,  a  wind  tunnel  test  program  was  conducted  on  an  inextensible, 
flexible,  semi-span  membrane  wing  in  the  Princeton  University  4'  x  5'  subsonic  tunnel.  The  model 
wing  was  constructed  from  two  thousands  thick  stainless  steel  and  was  of  triangular  planform  with 
an  aspect  ratio  of  9.8.  The  Reynolds  number  for  the  tests  was  6  x  10^.  The  model  was  designed  to 
allow  a  range  of  camber  and  washout  to  be  set  into  the  wing.  Observations  and  aerodynamic  data 
were  obtained  for  several  configurations  of  camber  and  washout.  This  report  contains  a  discussion 
of  the  wind  tunnel  tests  and  an  interpretation  of  the  results. 
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INTRODUCTION 

The  purpose  of  this  report  is  to  present  and  interpret  the  results  of  a  series  of  wind  tunnel  tests 
on  a  three  dimensional  inextensible  flexible  membrane  wing.  This  work  was  conducted  as  part  of  a 
basic  research  program  on  membrane  lifting  surface  aerodynamics  and  also  as  a  feasibility  study  to 
determine/evaluate  the  potential  effectiveness  of  the  wing  for  naval  use.  This  work  is  a  continuation 
of  earlier  studies  conducted  on  membrane  airfoils,  references  1  and  2. 

SCOPE 

The  report  will  discuss  the  results  of  a  series  of  wind  tunnel  tests  conducted  on  a  semi  span 
membrane  wing  mounted  in  Princeton  University's  4'  x  5’  subsonic  tunnel  and  their  interpretation. 
The  wmg  was  limited  to  a  triangular  planform  attached  to  a  steel  leading  edge  beam.  The  wing  was 
adjustable  to  produce  several  configurations  of  camber  and  washout  which  were  tested  over  a  wide 
range  of  angle  of  attack.  The  aerodynamic  data  obtained  consisted  of  lift,  drag,  and  leading  edge 
pitching  moments.  Qualitative  observations  were  made  on  unsteady  effects  at  both  low  and  high 
angles  of  attack. 


THE  WIND  TUNNEL  AND  WIND  TUNNEL  MODEL 

The  wind  tunnel  used  in  the  research  was  the  Princeton  University  4'  x  5'  single  return,  closed 
test  section,  subsonic  facility.  The  wind  tunnel  balance  uses  a  pneumatic/hydraulic  feedback  system 
to  measure  the  model  forces  and  moments.  The  wind  tunnel  model  consisted  of  a  semi  span  wing  of 
triangular  planform  with  an  aspect  ratio  of  9.8  made  of  sheet  stainless  steel,  two  thousandths  of  an 
inch  thick,  attached  to  a  leading  edge  spar  of  an  airfoil  sectioned  steel  tube,  see  figure  No.  1 .  The 
root  of  the  leading  edge  spar  was  mounted  in  a  spherical  bearing  and  the  tip  of  the  spar  was  attached 
to  an  adjustable  rod  support  off  the  main  beam  of  the  wind  tunnel  balance.  The  wing  root  trailing 
edge  attachment  was  connected  to  a  rod-end  spherical  bearing  mounted  on  a  manually  operated, 
screw  actuated,  chord  length  adjustment  slide. 

A  sketch  of  the  model  mounted  in  the  tunnel  test  section  is  shown  in  figure  No.  2  and  in  the 
photograph  figure  No.  3.  The  semi-span  wing  is  mounted  on  a  vertical  turntable  that  is  set  into  a 
reflection  plane.  The  turntable  could  be  made  to  rotate  remotely  to  adjust  the  angle  of  attack  of 
the  wing.  The  wmg  twist  and  camber  were  set  by  adjustments  made  at  the  wing  tip  and  at  the  wmg 
root  trailing  edge  location.  The  reflection  plane  was  mounted  on  a  r;gid  steel  bracket,  and  secured 
through  an  opening  in  the  wind  tunnel  floor  to  the  main  beam  of  the  wind  tunnel  balance. 

A  remotely  controlled  blockage  vane  was  installed  in  the  test  section  in  the  gap  between  the  re¬ 
flecting  plane  and  che  test  section  back  wall.  This  was  used  to  control  the  side  loading  on  the  reflec¬ 
tion  plane. 

HOW  THE  WIND  TUNNEL  TESTS  WERE  CONDUCTED 

The  wind  tunnel  balance  was  designed  to  operate  with  a  symmetrical  wind  tunnel  model  sup¬ 
ported  on  three  vertical  struts  adjusted  to  set  the  mode!  on  the  wind  tunnel  centerline.  In  the 
present  tests  to  increase  the  Reynolds  number  to  full  scale,  a  semi  span  membrane  wing  model  was 
built  mounted  on  a  large  reflecting  plane,  supported  off  the  main  beam  of  the  balance.  This  pro¬ 
duced  a  highly  non-symmetrical  loading  on  the  balance  system.  To  return  the  balance  system  to  an 
operable  range,  preloading  weights  were  added  to  several  of  the  active  balance  members. 
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After  balancing  the  model/balance  system,  the  system  was  calibrated  using  dead  weights. 

Preliminary  wind  tunnel  runs  were  made  to  measure  the  effects  of  the  reflection  plane  on  the 
balance  readings.  Using  a  dummy  wing  mounted  from  the  wind  tunnel  test  section  ceiling,  the  ef¬ 
fects  of  the  lifting  model  wing  on  the  tare  loads  of  the  reflecting  plane  were  measured.  These  results 
together  with  the  basic  tare  loads  for  the  reflecting  plane  were  used  to  determine  the  aerodynamic 
loads  on  the  wing  alone  from  the  wing  wind  tunnel  test  data. 

To  run  the  experiments,  the  membrane  wing  was  mounted  in  the  tunnel  and  the  root  excess 
length  adjusted  to  a  predetermined  setting.  The  root  excess  length  is  the  difference  between  the 
membrane  length  and  the  chord  length  at  the  root,  non-dimensionalized  by  the  root  chord  length 
and  expressed  as  a  percentage.  A  coarse  adjustment  of  the  wing  trailing  edge  shape  was  made  by 
moving  the  wing  tip  spar  mounting  location.  With  the  wind  tunnel  running  at  the  operating  speed 
q  =  10.5  LB.  FT*  and  the  wing  set  to  a  moderate  angle  of  attack  of  8:  to  10' ,  the  final  shape  of  the 
wing  was  set  by  adjusting  the  tip  support  -xternal  to  the  test  section.  A  major  controllable  variable 
having  a  significant  effect  on  the  shape  of  the  membrane  is  the  position  of  the  “Virtual  Apex."  Since 
the  wing  is  constructed  from  a  flat,  inextens;ble  but  flexible  membrane,  one  of  the  principal  radii  of 
curvature  at  any  point  on  the  surface  must  be  infinite.  This  condition  means  that  the  wing  is  made  up 
of  any  infinite  number  of  straight  lines  or  generator  lines  and  therefore  the  shape  must  be  part  of  a 
conical  or  cylindrical  surface.  It  should  be  noted  that  this  latter  constraint  is  not  as  restrictive  as  it 
might  seem.  The  basic  shape  of  the  cone  or  cylinder  does  not  have  to  be  circular  but  could  take  on 
shapes  that  include  reversal  of  curvature  so  the  membrane  still  has  the  freedom  to  assume  an  infinite 
variety  of  shapes.  The  definition  of  the  virtual  apex  is  given  in  figure  No  4  (a,  b.  c,  d  and  e).  Figure 
No.  4a  shows  the  general  definition.  Figure  No.  4  ib,  c.  d  and  e)  shows  large  and  small  virtual  apex 
lengths  and  infinite  and  negative  virtual  apex  settings. 

In  the  wind  tunnel  experiments  the  virtual  apex  lenqth  was  set  by  placing  a  telescope  on  a  line 
extending  out  from  the  leading  edge  of  the  membrane  and  viewing  a  zebra  striped  pattern  painted 
on  the  root  section  wall,  as  shown  in  figure  No.  3.  The  position  of  the  virtual  apex  was  set  by  adjust¬ 
ing  the  wing  tip  through  the  telescope,  i.e.,  the  virtual  apex  length  was  set  when  the  wing  was  seen  in 
edge  view  when  viewed  from  a  position  on  the  line  of  sight  of  the  membrane  leading  edge.  Because 
the  Reynolds  number  was  6x10®  and  close  to  the  boundary  layer  transition  number  of  5  x  10®  test 
runs  were  made  with  and  without  a  grit  boundary  layer  trip  strip  on  the  leading  edge  spar  with  the 
trailing  edge  tufted.  There  were  no  measurable  differences  in  any  of  the  forces  observed  nor  any 
visible  changes  in  the  flow  field.  It  was  concluded  that  the  critical  Reynolds  number  was  exceeded 
and  that  the  boundary  layer  was  essentially  fully  turbulent. 


3 


N  ADC-83 130-60 


After  setting  the  virtual  apex  position  and  a  given  root  excess  length  -  a  set  camber  -  lift, 
drag,  and  pitching  moment  data  were  recorded  along  with  a  pitot  tube  reading  of  the  tunnel 
dynamic  pressure.  Data  were  taken  for  a  range  of  angles  of  attack. 

The  inherent  wing  twist  effects  due  to  the  flexible  trailing  edge  of  the  wing  causes  the  wing 
surface  to  experience  areas  of  both  positive  and  negative  angles  of  attack  over  most  of  the  operating 
range.  This  point  is  discussed  in  detail  in  Appendix  A. 

In  preliminary  trial  runs  these  twist  effects  produced  some  unsteady  flutter  phenomena  in  the 
low  angle  of  attack  condition,  i.e.,  a  <  5°.  With  the  limitations  of  time,  funds,  and  scope  of  the  work 
a  decision  was  taken  to  limit  the  tests  to  angles  of  attack  a  >  5°.  In  unrecorded  test  runs  it  was  ob¬ 
served  that  the  wing  would  operate  down  to  a  =  -2°  for  many  of  the  conditions  studied  producing 
positive  lift.  The  only  exceptions  to  this  were  those  wings  with  large  root  excess  lengths. 

The  root  airfoil  shapes  observed  at  angles  of  attack  on  either  side  close  to  and  through  zero  were 
doubly  curved  ('S'  shaped)  and  closely  resemble  the  airfoil  shapes  seen  in  the  two  dimensional  work 
at  negative  angles  of  attack,  as  presented  in  Reference  1 . 

Unsteady  effects  were  also  observed  at  high  angles  of  attack  as  root  separation  occurred.  This 
had  little  effect  on  the  lift  of  the  wing  because  the  vibrations  were  held  to  small  amplitude  by  the 
active  lifting  outer  wing  span  which  acted  as  a  damper.  At  full  stall  of  the  wing  (i.e.,  as  all  sections 
stalled)  the  membrane  vibrated  in  a  controlled  way  at  a  frequency  dependent  on  the  root  excess 
length,  the  vertical  apex  length,  and  the  tunnel  dynamic  pressure. 

Table  No.  I,  page  28,  is  a  summary  of  the  wind  tunnel  test  program. 
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RESULTS 

The  results  of  the  wind  tunnel  tests  are  presented  in  figures  No.  5  through  20. 

Figure  No.  5  is  a  plot  of  lift  coefficient,  CL,  versus  angle  of  attack,  a,  for  several  experimental 
runs.  The  root  excess  length  and  the  virtual  apex  were  varied.  Runs  4,  7  and  9  have  a  similar  virtual 
apex  point  with  the  excess  length  varied.  Run  8  has  a  similar  excess  length  to  run  7  but  with  an  in¬ 
creased  virtual  apex  length. 

Generally  the  wing  configurations  with  the  larger  excess  lengths  stalled  at  higher  angles  of 
attack  and  produced  higher, maximum  lift  coefficients.  The  influence  of  the  excess  length  on 
the  maximum  CL  value  is  affected  by  the  position  of  the  virtual  apex.  This  can  be  seen  by  compar¬ 
ing  runs  7,  8  and  0.  Run  7  has  the  same  root  excess  length  of  2.899  percent  as  run  8,  but  run  7  has  a 
virtual  apex  of  12  FT  compared  to  18  FT  for  run  8.  Run  7  stalled  at  a  lower  angle  of  attack  and 
produced  a  smaller  maximum  lift  coefficient  than  run  8.  Run  0  with  an  excess  length  of  3.278  per¬ 
cent  and  a  virtual  apex  of  8  FT  produced  a  similar  lift  curve  to  run  7.  Decreasing  the  virtual  apex 
reduces  the  effect  of  the  root  excess  length  in  producing  lift.  Both  root  excess  length  and  the  posi¬ 
tion  of  the  virtual  apex  effect  the  twist  distribution  of  the  wing.  This  effect  is  discussed  in  ap¬ 
pendix  A.  The  effects  discussed  in  the  previous  comparison  of  runs  0,  7  and  8  are  due  to  the  effects 
of  twist  produced  by  the  various  settings  of  root  excess  length  and  the  virtual  apex  length.  (This 
effect  will  be  discussed  in  greater  depth  later  in  this  report.) 

It  can  be  seen  from  figure  No.  5  that  the  lift  curves  are  essentially  linear  with  the  slope, 

3CL/3ct,  decreasing  as  the  root  excess  length  decreases.  The  larger  excess  lengths  used  in  runs  7,  8 
and  0  of  2.899  percent,  2.899  percent  and  3.273  percent  respectively  produce  the  largest  lift  curve 
slopes.  The  lift  curve  slope  decreased  as  the  root  excess  length  was  reduced,  whereas  the  lift  pro¬ 
duced  at  5°  angle  of  attack  increased  as  the  root  excess  length  was  reduced.  The  only  exception  oc¬ 
curred  in  run  0  where  the  virtual  apex  was  8  FT  which  was  smaller  than  in  all  the  other  runs.  The 
short  virtual  apex  length  had  a  dominant  affect  over  the  root  excess  length  in  determining  the  wing 
twist  distribution  which  in  turn  affects  the  load  distribution  on  the  wing. 

It  should  also  be  noticed  that  the  lift  curves  for  run  0  and  run  8  when  extrapolated  back  to 
zero  lift  would  suggest  that  it  would  be  impossible  to  reduce  the  angle  of  attack  to  zero.  It  is  clear 
that  this  phenomena  associated  with  the  higher  excess  length  and  hence  higher  wing  twist  is  the 
reason  for  the  unstable  behavior  of  the  wing  at  the  small  angles  of  attack. 

The  results  of  the  lift  data  plotted  in  figure  No.  5  are  summarized  in  Table  1 1  on  page  28. 

Figure  No.  6  is  a  replot  of  data  runs  7,  4  and  9,  where  the  virtual  apex  length  is  held  approxi¬ 
mately  constant  at  12  FT.  The  plot  shows  the  effects  of  varying  the  root  excess  length,  XL,  on  the 
lift  characteristics  of  the  wing  independent  of  the  effects  of  the  virtual  apex  length.  From  Table  1 1 
it  can  be  seen  that  the  maximum  lift  coefficient,  CL,  increases  with  increasing  root  excess  length,  XL. 
The  angle  of  attack  range  increases  with  XL,  but  the  lift  developed  at  5°  angle  of  attack  is  lower  for 
increased  root  excess  length.  Figure  No.  7  is  a  replot  of  runs  7  and  8,  In  these  two  runs  the  root  ex¬ 
cess  length,  X  L,  is  held  constant  at  2.899  percent  and  the  virtual  apex  length  changed  from  1 2  FT  fcr 
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Run  7  to  18  FT  for  run  8.  For  the  same  root  excess  lengths  the  smaller  virtual  apex  length  produces 
smaller  twist  distribution  angles.  The  effect  of  the  smaller  twist  distribution  angles  is  to  lower  the 
maximum  lift  slightly  from  1.57  for  run  8  to  1 ,54  for  run  7  and  reduce  the  angle  of  attack  where  the 
maximum  lift  occurs  from  20°  for  run  8  to  18°  for  run  7.  The  slope  of  the  lift  curve  is  0.908/deg 
for  run  8  and  a  similar  0.810/deg  for  run  7.  The  major  difference  in  the  two  curves  is  that  the  curve 
of  run  7  is  displaced  vertically  above  the  lift  curve  of  run  8  This  results  in  a  coefficient  of  lift,  CL, 
of  0.535  at  a  =  5°  for  run  7  and  a  smaller  CL  of  0.370  for  run  8.  Since  the  root  excess  length  is  held 
constant  for  both  runs  7  and  8  and  the  virtual  apex  length  changed,  the  differences  in  lift  character¬ 
istics  between  run  7  and  8  are  due  solely  to  the  virtual  apex  length  change.  This  latter  result  is  in  con¬ 
trast  to  the  results  of  figure  No.  b  where  changing  the  root  excess  length  and  holding  the  virtual  apex 
length  produced  significant  changes  in  both  the  maximum  lift  attainable  and  in  the  lift  curve  slopes. 

The  major  effects  of  changing  the  root  excess  length  and  altering  the  virtual  apex  position  are  to 
change  the  spanwise  camber  and  twist  distribution  By  small  adjustment  of  the  root  excess  length  and 
the  virtual  apex  length,  various  mixes  of  camber,  and  both  geometric  and  aerodynamic  washout  are 
obtained.  The  effect  of  these  changes  on  the  wings  shapes  and  lift  characteristics  are  as  follows. 

Figure  No.  8  is  a  plot  of  the  theoretical  twist  angle  distributions  for  the  geometries  used 
in  runs  9,  4,  7  and  8  with  an  assumed  parabolic  camber  distribution  at  the  root.  Increasing 
the  root  excess  length,  holding  the  virtual  apex  constant  (runs  9,  4,  and  7),  increases  the  magni¬ 
tude  of  the  wing  twist  distribution  values  without  changing  the  shape  or  form  of  the  twist  distribu¬ 
tion  curve.  It  also  has  the  effect  of  increasing  the  camber  at  each  section  by  reducing  the  curvature  at 
each  section.  By  increasing  the  root  excess  length,  the  cone  that  the  wing  is  part  of  has  a  smaller  base 
circle.  The  aerodynamic  effects  of  increased  twist  and  camber  are  two  fold.  On  the  one  hand  the  in¬ 
creased  twist  allows  the  wing  to  operate  up  to  larger  angles  of  attack  and  produce  higher  maximum 
lift  values  before  tip  stall  occurs.  On  the  other  hand  the  increased  twist  produces  negative  angles  of 
attack  over  large  portions  of  the  wingspan  for  angles  of  attack  below  10°.  Although  at  the  larger 
angles  of  attack  there  was  evidence  of  stall  at  the  root,  the  lift  produced  by  the  tip  more  than  com¬ 
pensated  for  the  reduced  lift  at  the  root  due  to  root  Stall.  For  example,  from  the  wind  tunnel  test 
data  for  run  9  (XL  =  0.709  percent)  the  tip  stalled  at  1 1.5°  for  a  CL  MAX  of  1.131.  For  run  7  (XL  = 
2.899  percent)  the  tip  stalled  at  19°  for  a  CL  MAX  of  1.541. 

The  horizontal  line  drawn  across  the  graph  of  figure  No.  8  at  <t>  *  5°  indicates  the  areas  on  the 
wing  operating  at  positive  and  negative  angles  of  attack  for  a  root  angle  of  attack  of  a  =  5°. 

For  run  8  the  outer  76%  of  the  wing,  corresponding  to  61%  of  the  wing  area,  will  be  at  a  local 
negative  angle  of  attack.  Runs  4,  7  and  9  have  local  negative  angles  of  attack  for  the  outer  spans  of 
58%,  72%  and  46%  respectively,  corresponding  to  wing  areas  of  38%,  55%  and  25%  respectively. 

From  figure  No.  5  it  can  be  seen  that  the  reduction  in  the  CL  values  at  a  =  5°  for  runs  9,  4, 
and  7  is  an  example  of  how  the  increases  in  XL  with  the  accompanying  increases  in  twist  distribution 
results  in  the  outer  wingspan  sections  operating  at  negative  angles  of  attack  for  the  lower  angles  of 
attack  of  the  root  wing  section.  Note:  The  angle  of  attack  of  the  wing  is  defined  as  the  angle  of 
attack  of  the  root  wing  section. 

Comparing  runs  7  and  8  on  figure  No.  8  shows  that  increasing  the  virtual  apex  length  holding 
the  root  excess  length  constant  changes  the  shape  of  the  twist  distribution  curve  to  a  more  linear 
form.  These  changes  in  the  curve  produce  small  increases  in  the  twist  along  the  wingspan  without 
appreciable  changes  in  the  magnitude  of  the  twist. 

These  small  increases  in  the  wing  twist  values  produce  a  slight  increase  in  the  stall  angle  of 
attack  from  19°  to  20°  and  a  small  increase  in  the  maximum  lift  coefficient  from  1.54  to  1.57. 
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Decreasing  the  virtual  apex  length  produces  a  shorter  cone  which  increases  the  curvature  on  the  span- 
wise  wing  sections.  This  increase  in  curvature  has  the  effect  of  increasing  the  camber  on  every  wing 
section  excluding  the  root  and  tip  and  it  is  this  increase  in  wing  camber  that  accounts  for  the  vertical 
and  essentially  parallel  separation  of  the  curves  for  runs  7  and  8. 

Figure  No.  1 0  is  a  multiple  plot  of  coefficient  of  drag,  CD,  versus  angle  of  attack,  a,  for  the  ex¬ 
perimental  runs  0,  4,  7,  8  and  9  Comparing  the  drag  values  at  a  ~  5°,  the  smallest  drag  coefficient  is 
0.14  for  run  9  which  has  the  smallest  root  excess  length.  In  general  the  5°  drag  values  increase  with 
increasing  root  excess  length.  The  curves  for  each  run  are  of  parabolic  form  with  the  larger  root 
excess  lengths  having  a  broader  sweep  to  the  curve.  The  only  exception  to  this  rule  is  run  0. 

Run  0  has  the  largest  root  excess  length  of  3.273  percent  but  the  smallest  virtual  apex  length  of 
8  FT  (L  =  2.67).  From  figure  No.  9  it  can  be  seen  that  the  twist  distribution  curve  for  L  =  2.67  (a 
virtual  apex  of  8  FT)  is  lower  than  the  curve  for  L  =  6.0  (a  virtual  apex  of  18  FT)  so  even  though  the 
larger  excess  length  should  have  produced  the  largest  camber  and  twist  distribution  values,  these  ef¬ 
fects  were  completely  overpowered  by  the  effects  of  the  smaller  virtual  apex  length. 

A  comparison  between  runs  7  and  8  on  figure  No.  10  shows  the  effect  on  the  drag  of  varying 
the  virtual  apex  position  holding  the  root  excess  length  constant.  Run  8  with  the  longer  virtual  apex 
length  of  18  FT  (L  =6)  generates  less  drag  for  a  given  angle  of  attack  than  run  7  with  a  virtual  apex 
length  of  12  FT  (L  =  4).  This  indicates  that  as  the  wing  shape  becomes  more  cylindrical  it  produces 
less  drag. 

Figure  No.  1 2  is  a  multiple  plot  of  the  experimental  lift  to  drag  ratios.  All  the  curves  have  the 
same  general  form.  They  rise  for  increasing  angle  of  attack  to  a  peak  value  then  fall  off  to  the  sta'I 
value.  In  general  the  longer  root  excess  iength  gives  a  larger  L/D  ratio  and  this  occurs  at  a  larger  angle 
of  attack.  This  is  consistent  with  the  previous  analysis  of  the  lift  and  drag  curves  as  is  the  behavior 
of  run  0  which  does  not  fit  the  general  pattern.  It  might  be  expected  that  L/D  for  run  0  with  the 
highest  root  excess  length  would  peak  at  the  highest  angle  of  attack.  But  the  small  virtual  apex  length 
of  run  0  dominated  the  twist  distribution,  overshadowing  the  effects  of  increased  root  excess  length. 
A  comparison  of  runs  8  and  7  in  figure  No.  1 2  shows  that  the  iarger  apex  length  of  18  FT  ( L  =  6) 
for  run  8  versus  1 2  FT  ( L  =  4)  for  run  7  produced  an  increase  in  the  L/D  ratio  of  from  6  3  to  6.6 
and  these  values  occurred  at  a  =  16°  and  20v  respectively.  Hence,  the  larger  the  vertical  apex  length, 
the  more  cylindrical  in  shape  the  wing,  the  higher  the  maximum  lift  to  drag  ratio  obtainable.  Com¬ 
paring  the  L  D  values  at  a  -  5"  shows  that  L/D  ratios  decrease  as  the  amount  of  twist  increases. 

Referring  to  run  9  of  figure  No.  1 2  an  L/D  of  5.4  at  a  =  5'  is  obtained  with  the  smallest  root 
excess  length  ratio.  Whereas  runs  4,  7,  0  and  8  have  L/D  ratios  at  a  =  5"  of  4,  3,  2.6,  and  2.2 
respectively  showing  that  L,  D  decreases  as  the  twist  increases.  Next  refer  to  figure  No  9  where  the 
theoretical  twist  angle  distributions  are  plotted  for  runs  8  and  0.  Run  0  with  a  greater  excess  length 
than  run  8  has  a  smaller  virtual  apex  length.  This  reduces  the  magnitude  of  the  twist  over  90  percent 
of  the  semi  wingspan.  It  is  for  this  reason  that  run  0  with  the  largest  root  excess  length  does  not 
follow  the  pattern,  that  of  increasing  L/D  for  increasing  XL.  Again  the  small  virtual  apex  dominates 
the  root  excess  length  in  effecting  the  twist  distribution, 

The  experimental  lift  data  plotted  in  figure  No.  13  was  taken  to  give  some  indication  of  the 
virtual  apex  on  the  lifting  characteristics  of  the  wing. 

The  data  of  run  2  was  taken  at  a  constant  root  excess  length  of  1 .429  percent  and  at  a  constant 
angle  of  attack  of  103.  Data  from  run  4  was  plotted  for  a  reference.  At  an  XL  =  1 .429  percent  3t 
a  -  10'  the  wing  with  the  smaller  apex  length  has  the  higher  lift  coefficient.  The  lift  coefficient 
drops  as  the  virtual  apex  is  increased.  The  lowest  value  of  lift  coefficient  was  recorded  at  the 
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negative  virtual  apex  of  -3.6  FT  (L  =  -1.2).  Referring  to  figure  No.  1 1  (a  plot  of  the  theoretical 
twist  distribution)  it  can  be  seen  that  apex  length  increases  the  amount  of  twist  with  the  negative 
value  of  L  producing  the  greatest  twist  over  the  largest  amount  of  semi  wingspan.  These  progessively 
larger  values  of  twist  for  the  larger  excess  lengths  caused  progressively  larger  portions  of  the  semi 
wingspan  to  operate  at  negative  angles  of  attack  for  a  root  section  angle  of  attack  of  10°.  These 
results  are  consistent  with  results  plotted  in  figure  No.  5  where  increases  in  the  wing  twist  distribu¬ 
tion  reduce  the  lift  coefficient  for  angle  of  attack  below  12°.  At  low  root  angles  of  attack  the  twist 
on  washout  produced  locally  negative  angles  of  attack  for  large  portions  of  the  outer  wingspan, 
reducing  the  overall  lifting  capabilities  of  the  wing.  At  the  high  root  angles  of  attack  the  twist  al¬ 
lowed  the  outer  portions  of  the  wing  to  produce  unstalled  lift  that  improved  the  overall  lifting 
capabilities  of  the  wing. 

Figures  No.  14,  15,  16,  17  and  18  are  flow  visualization  studies  that  show  the  effect  of  wash¬ 
out  on  the  flow  over  the  wing  at  different  angles  of  attack.  The  wing  had  a  root  excess  length  cf 
1.429  percent  which  produced  a  camber  of  7.3  percent  of  the  root  section  and  had  a  virtual  apex 
of  1 1  FT  (L  =  3.67).  In  figure  No.  14  the  wing  is  at  an  angle  of  attack  of  14°  and  the  flow  is  fully 
attached.  In  figure  No.  1 5  the  wing  is  at  an  angle  of  attack  of  16°  and  there  is  intermittent  separa¬ 
tion  at  the  root.  In  figure  No.  16  a  =  20°  and  the  flow  is  separated  over  the  first  30  percent  of  the 
wingspan.  In  figure  No.  1 7,  a  »  24°  and  the  flow  is  separated  over  70  percent  of  the  wing.  In  figure 
No.  18,  a  =  26°  and  90  percent  of  the  wing  is  separated.  In  figures  No.  15,  16  and  17  the  tufts  are 
not  completely  parallel  to  the  stream  flow  direction  but  show  some  convergence  towards  the  reflec¬ 
tion  plane.  The  tufts  indicate  that  at  the  higher  angles  of  attack  the  lifting  wing  tip  feeds  an  inflow 
of  air  into  the  root,  which  tends  to  reduce  the  severity  of  the  wing  root  stall  and  to  allow  the  high 
cambered  airfoils  to  produce  substantial  lift  at  the  high  angles  of  attack. 

Figure  No.  19  is  a  crossplot  of  lift  coefficient  versus  root  excess  length  with  root  section  angle 
of  attack  as  a  parameter.  For  a  <  12°  the  curves  show  that  increasing  the  root  excess  length  reduces 
the  lift  coefficient  of  the  wing.  This  is  another  manifestation  of  the  increased  washout  produced  by 
increasing  the  root  excess  length.  Above  a  =  12°  the  washout  sections  of  the  wing  operate  at  a  posi¬ 
tive  angle  of  attack  so  that  for  a  >  12°,  increasing  the  root  excess  length  will  improve  the  lift 
coefficient. 

The  experimental  pitching  moment  data  is  presented  in  figure  No.  20.  The  form  of  the  curves 
are  similar  for  every  run  producing  a  consistent  family  of  curves.  The  pitching  moment,  CM,  is  neg 
ative  and  becomes  more  negative  in  all  cases  as  the  wing  root  section  angle  of  attack  increases.  The 
rate  of  decrease  is  not  linear  but  falls  off  as  a  increases  and  is  similar  for  all  the  runs.  Again  the  twist 
distribution  is  the  major  factor  in  analyzing  the  data.  In  general,  the  larger  the  root  excess  length 
the  larger  the  pitch  down  moment.  This  is  due  to  the  larger  inherent  pitch  down  couple  found  at 
zero  lift  for  a  cambered  airfoil.  The  larger  the  camber  the  larger  the  pitch  down  couple.  The  excep¬ 
tion  to  this  rule  is  found  in  run  8  in  which  the  virtual  apex  is  larger  than  in  all  the  other  cases.  Wnen 
run  8  is  compared  to  run  7,  at  the  same  root  excess  length,  the  wing  will  have  a  more  cylindrical 
shape  and  the  wing  twist  distribution  will  be  larger.  The  more  cylindrical  shape  will  produce  sec¬ 
tions  with  larger  curvatures  in  the  wing  which  means  smaller  cambers.  In  this  case  the  smaller 
camber  will  reduce  the  pitch  down  moments  at  the  lower  angle  of  attack.  But  the  larger  twist  distri¬ 
bution  produced  by  the  larger  virtual  apex  will  have  the  major  influence.  As  discussed  earlier  run 
no.  8  with  the  higher  twist  produces  lower  lift  at  the  lower  angles  of  attack  due  to  the  larger 
portion  of  the  outer  wingspan  washout.  It  is  this  lower  lift  at  the  smaller  angle  of  attack  that  contri¬ 
butes  to  the  low  pitch  down  moment.  It  can  also  be  seen  at  angles  of  attack  above  13°,  runs  7  and  8 
behave  identically.  This  is  because  at  the  higher  angle  of  attack  a  major  portion  of  the  lift  is 
produced  on  the  outer  wingspan. 
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CONCLUSIONS 

The  major  influences  in  determining  the  aerodynamic  characteristics  of  the  wing  shapes 
studied  were  the  root  excess  length  and  the  position  of  the  virtual  apex.  The  setting  of  the  root  ex¬ 
cess  lengtn  and  the  virtual  apex  controlled  the  root  camber  and  the  geometric  and  aerodynamic 
washout  throughout  the  wing. 

In  general,  the  higher  root  excess  lengths  produced  higher  maximum  lift  values  at  higher  angles 
of  attack.  A  reduction  in  the  virtual  apex  length  reduced  the  maximum  lift  values. 

The  maximum  lift  to  drag  ratios  for  all  the  configurations  studied  were  approximately  constant 
at  6.0  ;  0.5. 

Mildly  unsteady  behavior  was  observed  for  angles  of  attack  below  5'. 

For  a  given  root  camber  distribution  and  a  given  virtual  apex  position  the  wing  surface  is  com¬ 
pletely  defined. 

An  inextensible  flexible  material  can  be  used  to  produce  a  stowable  membrane  wing,  i.e.,  the 
wing  is  produced  from  an  initially  flat  surface  and  is  therefoie  capable  of  being  rolled  up. 

When  a  flexible  inextensible  membrane  is  formed  into  a  curved  surface  it  develops  relatively 
high  bending  stiffness  in  the  direction  perpendicular  to  the  curvature.  This  quality  could  be  used 
to  develop  a  wing  of  better  aerodynamic  shape  and  increased  area  since  a  wing  in  the  curved  state 
has  bending  stiffness  in  the  spanwise  direction. 

The  membrane  wing  model  has  demonstrated  its  potential  for  being  an  inexpensive,  stowable, 
easy  to  build  wing  that  could  have  many  possible  Navy  applications. 


RECOMMENDATIONS 

1.  Testing  of  other  shapes  including  planform,  aspect  ratio,  sweep  and  dihedral  variation  would 
lead  to  a  set  of  design  criteria  to  optimize  the  wing  design. 

2.  The  flutter  and  vibrations  that  appear  at  high  and  low  angles  of  attack  should  be  studied. 

3.  Tests  to  broaden  the  Reynolds  number  range  should  be  conducted.  Both  high  and  low  speed 
regimes  should  be  tested. 

4.  A  study  of  the  effects  of  camber  and  washout  in  the  overall  flowfield  around  the  wing  would 
be  of  great  practical  value  in  optimizing  the  aerodynamic  characteristics  of  the  wing.  The 
present  design  allows  camber  and  washout  to  be  changed  at  will  and  so  could  serve  as  a  model 
test  bed. 
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Figure  1.  Three-Dimensional  Wing  Semi  Span  Model 


re  2.  Wind  Tunnel  Model 
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e  3.  Model  in  Test  Section 
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(a)  GENERAL  CONFIGURATION 


LARGE  VIRTUAL  APEX 
(L  >10) 


(c)  SMALL  VIRTUAL  APEX 
(0<L<4) 

(L=  1  SHOWN) 
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(CYLINDRICAL) 
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(e)  NEGATIVE  VIRTUAL  APEX 

(L  <0) 


Figure  4.  The  Definition  of  the  "Virtual  Apex' 
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Figures.  Experimental  Lift  Curves 
Runs  0,  4,  7,  8,  9 
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Figure  7.  Experimental  Lift  Curves 
Excess  Length  Constant  Runs  7,  8 
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ANGLE  OF  ATTACK 


ANGLE  OF  ATTACK 
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xperimental  Lift/Drag  Cuives 
Runs  0,  4,  7,  8  and  9 


ANGLE  OF  ATTACK 
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Figure  18.  Flow  Visualisation  XL  =  1 .429  percent  (7.3  percent  Camber) 
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Figure  19.  Lift  Versus  XL  with  o  as  a  Parameter 
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ANGLE  OF  ATTACK 

Figure  20  Experimental  Pitching  Moment  Curves 
Runs  0,  4,  7,  8,  and  9 
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TABLE  I 


SUMMARY  OF  THE  WIND  TUNNEL  TESTS  RUNS 


— 

RUN  " 

r 

APEX 

XL 

a 

COMMENTS 

*0 

8  FT 

3.273% 

5—19 

1 

-3.6*“ 

1.429% 

10 

Vary  Apex 

!  2 

8  —  “  | 

1 .429% 

10 

Small  Apex 

3 

1  1  FT 

1.429% 

5  —  15 

Free  Boundary  Layer  Transition 

4 

1  1  FT 

1.429% 

5-16 

Fixed  Transition 

5 

Flow  Visualization  —  Photos  - 

-  Tufted  -  Fixed  Transition 

6 

Flow  Fisualizatio 

n  —  Photos  - 

-  Tufted  -  Free  T ransition 

7 

11  FT  | 

2.899% 

5-20 

Large  XL 

1  8 

18  FT  . 

2.899% 

5-22 

Large  XL,  Larger  APEX 

9 

12  FT 

1 

1 

0.709% 

5-12 

Small  XL 

TABLE  II 

SUMMARY  OF  AERODYNAMIC  LIFT  DATA 


RUN  = 

XL% 

ZR% 

FT 

APEX 

~ 1 
CL  @  a 

CL5° 

ifCL  do' 

9CL  da 

FLAT  PLATE 

8 

2  899 

10.7 

18 

wm 

•  . 

0.908  (2 jt) 

0 

3.273 

11.3 

8 

•  -  1 

0.783  (2n) 

7 

2.899 

10.7 

12 

H9 

0.535 

1 

0.810  <2w) 

4 

1.429 

7.3 

1 1 

1.328 

15s 

0.651 

0.073 

0.665  (2ji t) 

9 

0.709 

5.2 

12 

1.131 

11.5; 

0.750 

0.059 

0.535  (2 tt> 
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APPENDIX  A 


THE  TWIST  DISTRIBUTION  FUNCTION 


For  a  triangular  flat  surface  and  assuming  small 
angles 


p— — — 

x  =  x!  - 

-(£>  j 

© 

At  Point  P  on  the  trailing  edge. 
Az_  _  yv  —  v  _  1  -  y 
zR  yv  yv 

'  Az  =  zR  1 
1 _ 

(i 

V  yv  /  | 

© 

Also  Ax 

CR  ~ 

Yt  -  y 
yt 

Ax  =  CR 

i _ 

V  yt  ■ 

)|© 

m 


-60 
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From  equation  (?)  and  (?) 


zR 


Twist  Tan  o  =  Az  = 


('-£) 


Non  d:n’,e niMnalize  on  yt 


Define  '  =X;  —  =  Y;  —  =  Z;  =  N  ;  L  ~  — 

.  i  yt  yt  yt  yt 


(X)  Where 
VZR  =  f(XR) 


For  a  parabolic  airfoil  (Also  =  circular  arc) 


Leading  Edge 


•  • .  1 
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At  the  root  after  none  dimensionalizing. 


x,ix  I _ ! _ 

By  similar  triangles. 

XR  -  X  =  V  AND  X  =  X.  ( 1  — V)  from  equation  (T). 
Xfi  -  Xv  Yv  W 

XR  -  X  -  Y_  (XR  -  Xv) 

Yv 


ZR  =  4ZR  (  XR  -  Xo\ 

1  -  /  XR  -  Xo\ 

\  x1  / 

V  Xt  ) 

(  1  ~  Y\  = 

(  X  -  Xv  Y  ) 

V  Yv ' 

V  Yv  ' 

(x  -  Xv  y_\ 
' _ Yv  ■ 

/I  -  Y_\ 


Sut  rtituting  equations  (T)  and  (i^)  into  equation  (T)  gives  twist  angle  function  for  the  most 
general  case  of  a  parabolic  root  profile  with  a  linear  trailing  edge,  i.e.,  with  dihedral  and  sweep. 

For  the  Princeton  experiments  the  dihedral  and  the  sweep  were  zero.  With  these  restrictions 


Xo  =  Xv  =  0  ;  X1  =  N  And 
We  have 

ZR  =  4ZR  XR  /  1  -  XR\ 


•W 


X  =  N  ( 1— Y) 
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Equation  ("T)  becomes 


i  fY  (L  -  I)1 

=  42  R  r 

1  -  xrI 

L  L  -  Y  . 

N  L 

J 

But 


ZR  =  ?R  'cR  =  ZR  A 

~t—  — / —  -T7T-  b  ZRr  i.e,  maximum  camber  height  at  the  root,  non-dimensionalized 

M  yt,-  Vt  CR  u  ,  .  |  3 

on  the  chord  lengtn. 


Hence  0  is  not  a  function  of  wing  length  if  z  (the  camber  distribution)  is  non-dimensionali2ed  on 
the  chord  length  as  is  common  practice. 

I - 


TAN  o  =  42RC 


r y  (l  - 1) 
L  (L  -  Y) 


0 


Equation  is  the  twist  distribution  function  for  the  Princeton  triangular  airfoil  tests. 


A-5 
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APPENDIX  B 

HOW  THE  CAMBER  AND  VIRTUAL  APEX  LENGTH 
EFFECTS  THE  WING  TWIST  DISTRIBUTION 
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The  twist  angie  distribution  function  is  developed  in  appendix  A.  The  function  assumes  a 
parabolic  wing  root  profile,  a  reasonable  approximation  of  an  actual  wind  tunnel  test  profile.  This 
function  is 

<t>  =  TAN'1  4ZRCY  (L  -  1) 

<L-Y) 


where 


ZRC  is  the  maximum  profile  camber  non-dimensionalized  on  the  root  chord 
Y  is  the  spanwise  location  non-dimensionalized  on  the  wing  semi-span 
L  is  the  position  of  the  virtual  apex  non-dimensionalized  on  the  wing  semi-span 
0  is  the  angle  of  twist 


For  2r  <  0.1,  <p  is  linear  in  ZRC  to  within  1  percent.  Hence  the  twist  angle  distribution  is  essentially 
proportional  to  the  camber  for  any  value  of  L. 


Figure  No.  1 1  is  a  plot  of  angle  of  twist,  0,  3gainst  percentage  semi-wingspan  location,  y,  with 
the  virtual  apex  position,  L,  as  a  parameter.  As  L  approaches  infinity  (L  -*•  <»)  the  wing  surface  be¬ 
comes  part  of  a  cylinder  and  TAN  0  =  4ZRC  y.  In  this  case  the  twist  distribution  is  linear  in  both 
camber,  2RC,  and  wingspan  location,  y.  As  L  approaches  unity  { I_  — *•  1 )  the  wing  surface  becomes 
part  of  a  conical  surface  with  the  apex  of  the  cone  at  the  wing  tip  position.  In  this  case  both  the 
leading  and  trailing  edges  of  the  wing  are  co-planar  and  TAN  0  =  0  which  implies  that  there  is  no 
twist  distribution  in  the  wing.  In  practice  it  is  impossible  to  achieve  zero  twist  distribution  with  a 
flexible  trailing  edge. 


For  a  fixed  straight  leading  edge  spar  the  virtual  apex  length  cannot  occur  between  the  limits 
of  one  and  zero. 

For  a  large  negative  virtual  apex  position  the  wing  surface  again  becomes  part  of  cylinder  and 
as  in  the  case  of  large  positive  L  the  twist  distribution  is  linear  in  both  camber,  ZRC,  and  wingspan 
location,  y. 

For  negative  L  values  close  to  zero  the  twist  distribution  angle  approaches  a  constant  value, 
that  of  the  limiting  value  at  the  wing  tip,  A  value  which  is  completely  determined  for  a  given  root 
profile  shape.  For  the  selected  representative  example  of  a  10  percent  parabolic  profile  this  limited 
value  of  the  wing  tip  twist  angle  is  2 1 .80° . 


B-2 


